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Abstract 

A nonlinear kinetic chemotaxis model with internal dynamics incorporating signal transduction 
and adaptation is considered. This paper is concerned with: (i) the global solution for this model, 
and, (ii) its fast adaptation limit to Othmer-Dunbar-Alt type model. This limit gives some insight 
to the molecular origin of the chemotaxis behaviour. 

First, by using the Schauder hxed point theorem, the global existence of weak solution is proved 
based on detailed a priori estimates, under some quite general assumptions on the model and the 
initial data. However, the Schauder fixed point theorem does not provide uniqueness. Therefore, 
additional analysis is required to be developed to obtain uniqueness. 

Next, the fast adaptation limit of this model is derived by extracting a weak convergence 
subsequence in measure space. For this limit, the first difficulty is to show the concentration effect 
on the internal state. When the small parameter e, the adaptation time scale, goes to zero, we 
prove that the solution converges to a Dirac mass in the internal state variable. Another difficulty 
is the strong compactness argument on the chemical potential, which is essential for passing the 
nonlinear kinetic equation to the weak limit. 
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1 Introduction 

Chemotaxis is a mechanism by which cells or bacteria efficiently and rapidly respond to changes in 
the chemical composition of their environment, for example approaching chemically favorable environ¬ 
ments or avoiding unfavorable ones. This behavior is achieved by two major steps: (i) detection of the 
signal and (ii) integrate the signals received from receptors that triggers the response. We consider the 
kinetic chemotaxis model with internal dynamics incorporating signal transduction and adaptation 

dtp + v- VxP + = / T{v, v', m)p{x, v',m, t) - T{v', v, m)p{x, v, m, t)dv', 

Jv 

' —AS + S = n{x,t):= f f p{x,v,m,t)dvdm, 

Jo Jv 

p{x, v,m = 0,t) = 0, p{x, v,m = -|-oo, t) = 0, 
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where p{x,v,m,t) denotes the cell density at position x € velocity v & V, with internal state 
m > 0 (for example, methylation level) and time t, S{x,t) is the chemical potential. The velocity 
space V is assumed to be a compact domain in W^, typical example being unit ball {|u| < 1}. The 
parameter e in the model characterises the adaptation time scale. 

This model has been developed in [3ll [3l [27], and further elaborated by [El [E], to connect the 
aspects of the signal transduction and response to the macroscopic equations by using moment closure 
techniques. This type of model was also studied in some other references, for example [SaiMllSTI , to 
bridge the molecular level pathway dynamics with the cellular behavior such as cell population level 
motility in some biological systems. These works made possible the development of predictive agent- 
based models that include the intracellular signaling pathway dynamics. For instance, it is of great 
biological interest to understand the molecular origins of chemotactic behavior of E. coli by deriving 
a population-level model based on the underlying signaling pathway dynamics. We further mention 
that, this model also has close relation to the kinetic theory for active particles [6l[5| (KTAP), which 
has been applied to model various complex systems in life sciences, for instance biological growing 
tissues in [22] , social systems |7] , behavioural economy [T] or epidemics with gene mutations m- 

The mathematical study of chemotaxis dated back to as early as the work of Patlak [29] in 1953, 
and further Keller and Segel derived one of the best studied models in mathematical biology at the 
macroscopic population level [231 [Ml Eg. We refer m and references therein on this model. The 
famous Keller-Segel model has successfully explained chemotactic phenomena in slowing changing 
environments [36]. On the other hand, in order to understand bacterial behavior at the individual 
level, the Othmer-Dunbar-Alt kinetic model was proposed by Alt (1980) [2] and Othmer et al. (1988) 
[28] for the description of the chemotactic movement of cells in the presence of a chemical substance, 
and the Keller-Segel model can be derived by taking the hydrodynamic limit of kinetic models (see 
[El [El [El El] and references therein). In the conventional Othmer-Dunbar-Alt kinetic models, a 
(biased) velocity jump assumption in the turning kernel is always assumed. 

This paper is concerned with the study of the kinetic chemotaxis model with internal dynamics 
Our first goal is to prove existence of global solutions for the Cauchy problem of (jl.ip under 
some quite general assumptions. We assume the initial condition 

Po{x,v,m) >0, (1-1- {x))po{x,v,m) € T^(M^ x K x po{x,v,m) € L°°(Mf x K x (1.2) 

where {x) = vT+lE, and 

poo 

Po{x,v) := / pQ{x,v,m)dm £ X V). (1.3) 

Jo 

We will use the notation 

poo 

p{x,v,t) := / p{x,v,m,t)dm. 

Jo 

The assumptions on F{m,S) and T{v,v',m) are: 

• Assumption on the adaptation rate F. There exists a non-negative, non-decreasing continuous 
function n(.) G C(M) such that 

|5^E(m,5)| <n(5). (1.4) 

• Assumption on the turning kernel T. The turning kernel T is positive and uniformly bounded: 
there exist a positive constant C-p such that 

0 < T{v,v',m) < C-p. (1.5) 
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Under the above assumptions, we prove the global existence of weak solution to the system (jl.ip . see 
Theorem 12.21 in Section [2l However, the Schauder fixed point theorem does not provide uniqueness. 
With further assumptions on the coefficients and initial data, the regularity and uniqueness of the 
global weak solution are also considered in Section [3l 


The second goal is to study the fast adaptation limit as e ^ 0 of dni). For the reasons below 
(compatible with |32t [33] ), we make further assumption on F that 


There exist positive constants m± and mo{S),m- < mo{S) < m+, such that 
F{m,S)>0 for m<mo{S), F{m,S)<0 for m>moiS). 

And to simplify the proof, we also assume the initial data does not concentrate on large m, i.e., 

POO P P 

/ / / mpodvdxdm < oo. 

Jo Jm Jv 


( 1 . 6 ) 


(1.7) 


Under the above assumptions, we study in Section (U as e —)• 0, how the kinetic chemotaxis model 
with internal dynamics (II.ip is close to the Othmer-Dunbar-Alt type kinetic model 

dtp + v-V:;cp = T{v,v',moiS))p{x,v',t)-T{v',v,mo{S))p{x,v,t)dv', 

Jv 

—AS + S = n{x,t) := / p{x,v,t)dv. 

Jv 

The resulting limit. Theorem 14.31 gives some insight to the molecular origins of the chemotaxis be¬ 
haviour. 


We conclude this introduction by mentioning that, other types of scaling and related limits can be 
considered. For example, based on both fast adaptation and stiff response, the paper m studied how 
the path-wise gradient of chemotactic signal arises from intra-cellular molecular content. See further 
remarks in the conclusion section. 


2 Global existence of solutions 


Existence of solutions for this type of problem has been considered in M- See also [9]. Here, we 
extend the results under more general assumptions. We also give uniform bounds which will be useful 
for later purpose. The parameter e does not play a role here thus we take e = 1 for simplicity and 
rewrite the original equation for p as 


dtp + V ■ VxP + dm[F{m, S)p] = / T{v,v',m)p{x,v',m,t) — T(v',v,m)p{x,v,m,t)dv', (2.1) 

Jv 

coupled with the steady elliptic equation for the chemical signal S 


— AS + S = n{x, t) 

:= / / p{x,v,m,t)dvdm, 

Jo Jv 

(2.2) 

and the boundary constraint 



p{x, v,m = 0,t) = 0, 

? 

II 

+ 

8 

II 

O 

(2.3) 
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Note that from Equation (|2.2I) . we have 


S{x, t) = G * n{x, t), 


where G is the Bessel potential. For later use, we recall some properties of the Bessel potential in 
below. 

Proposition 2.1 (Properties of Bessel potential 135^ ) 

(PI) G{x) > 0, ||G||ii(K.) = 1, 

(P2) For any spatial dimension d, 


G{x) G VI < g < 

(P3) For any a G M, /? > 0, 


(P4) For any spatial dimension d, 


d 


d-2' 


{we understand 


d 


d-2 


= oo when d = 1, 2), 


f |a;|"G^(x)hx < oo, 
J\x\>l 


VG(x) G L\R'^) 


We briefly mention that, for d = 1, G{x) = For d = 2, G{x) has the singularity of log\x\ at 

l^l —)• 0. In the case h > 3, the only singularity of G{x) at x = 0 behaves as 

G{x) ~ as \x\ —> 0. 

Also note that G{x) decreases exponentially as |x| —>■ oo, therefore we obtain (P2)-(P4). 

2.1 Statement of the result 

Theorem 2.2 (Global existence) Under the assumption (|1.4p on F, (II.5p on T, consider system 
dUD-dlS]) with initial data po satisfying (nai-dni); there exists a global weak solution {p, S) that 

\\Si;t)\\L^nL^iRi) < WpoWl^ + Vu\\po\\L^{l + 2VdGTte^^<^^rt^, Vt > 0, (2.4) 

lll'(') 'i ■) 0IILinL°°(K^xy xk;1i) — IIpoIIli + ||po||l°° ' ) > Vt > 0, (2-5) 

where Vd is the volume of V, C > 0 is a constant depends on the estimate of S, Cj- and initial data. 
Moreover, 

l|p(v,t)llLoo(K.xy) < \\po\\LAl + 2VdCTte^^'^^^^), Vt > 0, (2.6) 

l|n(-,t)||Loo(„.) < Frf||po||L-(l + 2VdCrte^^‘^^rt^, Vt > 0. (2.7) 

We recall that the global existence result was considered in CD, with more precise examples on T 
and F. They take the turning kernel T as a product of the turning frequency A and the kernel K, i.e., 

T{v,v,m) = X{m)K{v,v' ,m), 

in which the kernel K is non-negative and satishes the normalisation condition 

/ K{v, v', m)dv = 1, 

Jv 
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where ^ is a symmetric compact set in The simplest example of K is 

1 


K{v, v', m) = 




to assume constant turning probability of changing velocity from v' to u, or 


K{v,v',m) = k9, cos{6) = 


V ■ V 




to assume that the turing kernel is a function of the angle between original and new velocities. More 
general, one can assume that it is uniformly bounded by a constant , i.e., 

\K{v, v', m)\ < C. 

For the (output) turning frequency A(m), which is related to the (input) signal function seen by a cell, 
to prevent formation of singularities, the authors in m suppose the growth estimate 

X{m) <C\l + A{S) + \— + v-—\y \/x, V, t, 

where A S C'(M^) is a non-negative, non-decreasing continuous function. An assumption that, as we 
see it below, is not necessary for our mathematical treatment. 


2.2 A priori estimates 

• bounds on p, p and n. 

Take the integration of (j2.1h with respect to x, v, m variables, we have the mass conservation 


A 

dt 



p{x, V, m, t)dvdmdx = 0, 


which is also 


IIp(') ■) ■5^)lll,i(Kdxyx]R+) “ ■) ^)IIl1(K^xV") — lll,i (R^) 



IV 


Po{x,v, m)dvdmdx 


IIpoIIli. 


( 2 . 8 ) 


• L°° bounds on p and n. 

On the other hand, take the integration of (j2.ip only with respect to m over [0, oo) and use the 
boundary condition (j2.3p . we have 


poo p 

dtp + v-'VxP= / / T{v,v',m)p{x,v',m,t) — T{v',v,m)p{x,v,m,t)dv'dm, 

Jo Jv 

then it can be further represented by 


p{x,v,t) 


Po{x 


pt poo p 

vt,v)+ / / \^T{v,v',m)p{x 

Jo Jo Jv 

— T{v', V, m)p{x 


v{t 

v{t 


s),v', m, s) 

s),v, m, s)'^dv'dmds. 
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Recall the initial condition (II.3p and assumption (jl.Sp on T, take the x V) norm on both sides 

of the above equation to get 

IIM') < IIpo||l°° + 2VdC'r [ ||p(-, •, s)||2,oo(iRdxv)ds. 

Jo 

Then, apply Gronwall’s inequality, we readily find p2.6l) . 

Next, from the definition of re, we have that 

thus we have (12.71) . 

• and L°° bounds on S. 

Recall that S{x,t) = G * n{x,t), by Young’s inequality for convolution, it is direct that 

ll'S'(-,i)||Li(Rrf) < l|G'llLi(R^)ll^('!OllLi(Rrf) = IIPoIIli, (2-9) 

where (12.Sp is used, and 

ll‘S'(-T)llL°°(Rrf) < l|G'llLi(R^)lh('j^)llL°°(Rrf) < krf||po||L°°(l + 2VrfC'rte^^‘^‘^'^*), 
in which (j2.7n is used. Then this inequality and p2.9l) gives the estimate p2.4l) . 

• L°° bound on p. 

Now we introduce the characteristics of Equation (j2.ip by 


^ = v, ^ = 0, ^ = F(M(s),5(X(s),s)), (2.10) 

as as as 

and along time-backward characteristics starting at {x,v,m,t), we have for 0 < s < f, 

X(s) = X — u(f — s), M(s) = m — f T(M(r), S'(X(t), r))dT, 

J S 

integrate m along the characteristic from 0 to f we get 

p{x, V, m, t) = po(X(0), u, M(0)) - f dmF{M{ t), S(Kir), t))p{'K{t), v , M(t), T)dT 

Jo 

+ f f rpu, u', M(T))p(X(r), u', M(t), r) — T(u', u, M(T))p(X(r), u, M(r), T)(iu'(iT. (2-11) 
Jo Jv 

Now, take the x R x M+) norm on both sides, and under the assumptions (|1.4p - (ll.5l) . we have 




L°°(R^xyxR+) 


< IIpoII 


L°°(R^xyxR+) 


+ [n(maxJ5(X(s),s)|)+2RdCr] ^ IIM',T)llLoo(RdxyxR+)'^R 

applying Gronwall’s inequality again, we obtain (|2.5I) . 


( 2 . 12 ) 
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2.3 Proof of Theorem 12.21 

We use the Schauder fixed point theorem to prove global existence. 

Theorem 2.3 (The Schauder fixed point theorem ]ltA \31^ ) Let X be a normed vector space, and let 
K C X be a non-empty, bounded, and convex set. Then for any given continuous mapping f : K ^ K, 
with 4'{K) being pre-compact, there exists a fixed point x £ K such that 4>{x) = x. 

To apply this theorem, we separate the presentation into three steps. 

Step 1. Setting of the function space and the mapping. Fix Tq > 0, we define 

X:=LH[0,ro]xM^), 

and take a bounded convex subset K m. X defined by 

K = {S£ X\S > 0, \\S{;t)hi < IIpoIIli, l|5(-,t)||Loo < Prf||po||L-(l+2PdCrroe2^^^^^°), Vt e [0,ro]}. 

Then we dehne the mapping cf. We start from a function S £ K, and construct p £ L°°([0, Tq], x 

V X R(n)) according to (12.Ih with S fixed. This is possible because it is a linear operator and the 
characteristics are well defined according to ()2.10l) . 

Next, the integration of p with respect to {v,m) dehnes n £ L°°([0, Tq], L^(M^)). By the above a 
priori estimate (j2.8p and (j2.7l) . we have 

ll^^(•,^)llI,l(Kd) = IIpoIIli, vt £ [o,ro], 

l|n(-,t)||Loo(R.) < VfilPollL^(l + 2V,Crte^^‘‘^rt)^ Vf G [0,ro]. 

Note that by interpolation, for any q £ (l,oo), 

(2.13) 

Further, we construct S as the solution to 

poo p 

— AS + S = n := / / p{x,v,m,t)dvdm. (2-14) 

Jo Jv 

By using Besssel potential and Young’s inequality for convolution (as the proof of (12.9p ). we have 

l|S(-,t)||Li < IIpoIIli, l|S(-,t)||L- < Vfi\po\\L^{l + 2VdCTToe^^^^rTo^, Vf G [0,ro], 

thus it holds that T, £ K, therefore we have defined a mapping on K that 

(/> : 5 !-)• S. 

Step 2. Continuity. To prove the continuity of the mapping : S T,, we observe that, firstly, 
the mapping 

</>! : p I—)• S 

dehned by (I2.14p is continuous, because it is a bounded linear operator from to for 1 < q < oo. 
Secondly, the mapping 

(f)2 ■■ S p 
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is also continuous, which can be seen exactly from the representation formula (I2.1ip (and also (|2.12p i. 
Indeed, the characteristics are uniquely defined and continuous with respect to parameters, although 
they are not Lipschitz, neither relevant for DiPerna-Lions theory |12] . 

In conclusion, the composition (/> : S' i—>• S is continuous. 


Step 3. Compactness of the mapping. To use the Schauder fixed point theorem, we need further 
to prove that 

(p{K) is pre-compact in K. (2-15) 

The proof of (I2.I5P consists of the following three claims. 

Claim 1. Local compactness in space: 

vte[o,ro], vi < g < oo. 

This claim is clear by using standard elliptic regularity estimate |19] . since S is defined by (j2.14p and 
recall ()2.13l) that n{-,t) S L'^(]R^), for all q> 1. 

Claim 2. Local compactness in time: 


VtG[0,ro], Vl<q<oo. 
This is because, take time derivative of the elliptic equation for S we have 


-AdtU + dtS = dtn = 


roo r 

Jo Jv 


vpdvdm. 


Note that for any t, the integration on the right hand side of above is in L'?(M^), for any q > 1, the 
standard elliptic regularity theory shows that dtT, G W^’'^(R^), and actually 

\\dtT,{-,t)\\wi,<i(Rd)<C\\ / / up < C diam(y) II / pdv\\L^(^d^ 

Jo Jv Jv 

< C II / p(iu||^inLoo(]Rd) < C ||p||LinL°°(Kdxy)) ^ [O^^o]) 

J 

where C independent of S' is a constant given by elliptic estimate and C is a genetic constant which 
differs from line to line. 


Claim 3. Control for |x| ~ oo: 

/ {x)T,dx < oo, Vt G [0,To], where (x) := i/l + \x\‘^. 
jRi 

For this, note that S satishes (|2.14l) . then 

/ {x)A,dx = / {x)ndx + / {x)AT,dx, 

jRd JlRd JlRd 

in which the second term is bounded because 

[ (x) ■ AT,dx = [ A{x) ■ T,dx = [ r—h -ttt) ■ ^dx < dUS 

jRi JRi JRi \x) {xr 


(2.16) 


(2.17) 

















To control the first term on the right hand side of ()2.17l) . we multiply Equation (12.ip by (x) and 
integrate: 


A. 

dt 



iv 


{x)pdvdmdx + 



V ■ Vxpdvdmdx = 0. 


Then 


— [ {x)ndx = — [ [ [ {x)pdvdmdx = [ f f pdvdmdx < ||po||li- 

dt jRd dt jRd Jq Jv jRd Jo Jv \x} 

Therefore from the initial date dEl we have the bound 

/ {x)ndx < / {x)nodx+ t\\po\\Li = ||(x)po||li +i||po||Li, Vt G [OjTo], 

jRi JRi 

then we have proved (12.161) . and further, 

(x)S G X = L^([0,ro] X Mf). 


(2.18) 


The property (I2.18P yields: 

Veo > 0, 311 C [0,To] x bounded, measurable, such that ||E||j;^iqq j’pjxRd \n) < ^o- 
This is because 


r-To 


Tjdxdt < 




1 r 

W)Jo A 


\x\>R 


{x)T,dxdt —)• 0 as i? —)■ oo. 


(2.19) 


In conclusion, for any fixed 11 with finite measure defined in (12.191) . by Claim 1 and Claim 2, we 
have 

{S|q} is precompact in L^(H), 


together with (|2.19l) . the weighted in space control away from H, we proved (ri.lSp . by using strong 
compactness criterion in Brezis [8] (Theorem 4.26, Corollary 4.27, p. 111). 


Now we use the Schauder Theorem 12.31 to get a fixed point which proves existence of solution 
S G K C X = L^{[0,To] X R^) for all Tq > 0. And next, we can also construct the global solution 
p G n L°°(R^ xVx R^)) for all Tq > 0 by the above a priori estimates. The proof of 

Theorem 12.21 is complete. □ 


3 Regularity and uniqueness of the solution 

In above section, the existence of global weak solution is derived by using the Schauder fixed point the¬ 
orem. However, the Schauder fixed point theorem does not provide uniqueness. Therefore, additional 
analysis is required to be developed to obtain uniqueness. 

Note the property (Hd) of G in Proposition 12.11 and the a priori estimate (12.71) . by using Young’s 
inequality we have 

||9xi5'||L°o(iRd) < ||'92;.G|hi(Rd)||n||j;^oo(Rd) < oo. (3.1) 

Also recall (j2.4p that S{-,t) G L°°(R^), then we obtain 
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Proposition 3.1 Let S be the weak solution of (j2.ip - (l2.3p defined in Theorem \2.2l Then 

Si-,t) G 

By further assumptions on regularity of coefficients and initial data, we can get higher regularity of 
p. For this purpose, we further assume 

• There exists a non-negative, non-decreasing continuous function n(-) G C(M) such that 

\dsF{m, 5)1 -F \dmmF{m, 5)| -h \dmdsF{m, 5)| < n(5). (3.2) 

• There exist a positive constant Cj- such that 

|V„r(u,t>',m)| -I- \Vy'T{v,v',m)\ + \dmT{v, v', m)\ < Cr- (3.3) 


Proposition 3.2 Under the assumptions on Theorem \2.B . Further assume (13.2p - (13.31) and let initial 
data po G n VF^’°°(M^ x V x M+). Then the weak solution of (j2.ip - (j2.3p satisfies 


1,1 r>, TT/1,00/TiDfi 


p(-,-,-,t) G IT^’i nlF 


X P X 


(3.4) 


Proof. 1. Regularity on x. Differentiate (j2.1[) with respect to Xj, 1 < i < d, we have 

dtdxiP + V ■ VxdxiP + F{m, 5)5^5^,p = -2dmF{m, S)dxiP - dmdsF{m, S) dx^S p 

+ / T{v,v',m)dxiP{x,v',m,t) — T{v',v,m)dxiP{x,v,m,t)dv'. 

Jv 

Integrate along the characteristics defined by (j2.10p from 0 to t. 

dxiP{x,v,m,t) = (9^,po(X(0),u,M(0)) - f 29^F(M(r), 5(X(r), r))(9^,p(X(r), u, M(r), r)dr 

Jo 

- [ dmdsF{M{T),S{X{T),T))dxiS{X{T),T)p{'X{T),V,M{T),T)dT 

Jo 

+ [ [ r(i;,t>',M(r))9a,.p(X(r),i;',M(r),r) - T(u',u,M(r))9a;.p(X(r),t>,M(r),r)dt>'dr. 

Jo Jv 

Using the assumptions on the coefficients and initial data as stated above, we get 

\dxiP{x,v,m,t)\ < |9a;.po(X(0),t>,M(0))|-h 2n(max |5(X(s),s)|) f |9a,.p(X(r), u, M(r), r)|d7 

0<S<t Jq 

-hn(max |5(X(s),s)|)||9„,5||^cx,(]Kd) f lp(X(r),u,M(r),r)|dr 

0<s<t Jo 

FTVdCr [ |9xiP(X(r),u',M(r),r)|dr, 

Jo 


that is. 


\dxiP{x,v,mfi)\ < |5xiPo(X(0),u,M(0))| 

-h2[n( max |5(X(s),s)|)-F PdCr] / |9xiP(X(r), u, M(r), r)|d7 

0<s<t Jq 


(3.5) 


-hn(max |5(X(s),s)|)||9„,5||icx>(Rd) [ |p(X(r),u,M(r),r)ld7 

0<S<t ^ Jq 
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2. Regularity on v. Differentiate m with respect to r;*, 1 < i < d, we have 

+ V • + F{m, S)dmdy^p = -dmF{m, S)dy^p - d^iP 

+ / T{v,v',m)di;iP{x,v',m,t) — T{v',v,m)dy^p{x,v,m,t)dv' 


'V 


+ / dy^T{v,v',m)p{x,v',m,t) — dy^T{v',v,m)p{x,v,m,t)dv'. 

Jv 

Similarly, integrate along the characteristics defined by ()2.10p from 0 to t, and use the assumptions, 
we obtain 

\dy^p{x,v,m,t)\ < |9„,po(X(0),u,M(0))| + / |(9a;iP(X(T), u,M( t), r)|(iT 

Jo 

+n(max |S’(X(s),s)|) [ |(9^,p(X(r),u,M( t), r)|(ir 

0<s<t ./n 


+2VdCr [ \p{'X{T),v,M{T),T)\dT FVdCr [ |5„,p(X(r), u, M(r), r)|(i7 
Jo Jo 


3. Regularity on m. Next, differentiate (12.11) with respect to m, we have 

dtdmP “1“ X • ^xdmP “1“ F(j 71, S'jdmdmP — dmmF(fn, S')p 2dfyiF(j7l, S'jdmP 

+ / T{v,v',m)dmp{x,v',m,t) — T{v',v,m)dmP{x,v,m,t)dv' 

Jv 

+ / Tm{v,v',m)p{x,v',m,t) — Tmiv',v,m)p{x,v,m,t)dv'. 

Jv 

Integrate along the characteristics defined by (I2.10p from 0 to t, we obtain 
\dmP{x,v,m,t)\ < |(lmPo(X(0),U,M(0))| 

-[2n(max |5(X(s),s)|)+214C'r] [ |i9mP(X(T),u,M( t), r)|(i7 

0<s<i Jo 


(3.6) 


+ 


(3.7) 


+ 


[n(max |S’(X(s),s)|) + 21/rfCr] [ |p(X(t), u,M( t), r)|(iT. 

0<S<t Jn 


In summary, recall the estimate (12.5p for p and Proposition 13.11 for S, combining (|3.5l) - (l3.7p and 
taking L*? norms, using Gronwall’s inequality, we deduce, for 1 < g < oo, V t > 0, 

l|VxP(-,-,-,i)||L'i + ||V^p(-,-,-,t)||M + \\dmP{-,-V:t)\\L‘i < C{\\S{-,t)\\wr^,\\Po\\w^,<!,Vd,CT,Cr) < oo. 
Together with (12.5p . we have (13.41) . □ 


For later use, we derive the following a priori estimates. 

Proposition 3.3 Under the assumptions on Proposition \3.^ and further assume the initial data 
satisfy dmPo £ m)- Then for all t > 0, 

IIpIIl-(l1 + \\dmP\\L^(Lf^) < OO. 
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Proof. Firstly, recall the a priori estimate (12.6p on p thus p € L!^.y(L^). Note V is compact, then 
we readily have 

\\p\\l^{lI^) < oo. 

Next we prove pm £ take L’^{Ll^) norm on both sides of (13.71) . then the 

result follows by using again Gronwall’s inequality. □ 

We now state the uniqueness result. 

Proposition 3.4 Under the assumptions on Proposition \3.3\. the weak solution of ()2.ip - ()2.3l) is 
unique. 

Proof. We follow the steps in [26]. Assume both (pi, 5i) and {p 2 , S 2 ) be weak solutions of (I2.ip - (|2.3I) 
with same initial data pq satisfying (USD-Jop. Denote 

P = pi-P2, S = 81 - 82 . 


Then we have 


dtp + v- VxP + dm[F{m, 8 i)p + (F(m, Si) - F(m, 82 ))p 2 ] 



= / T{v,v',m)p{x,v',m,t) — T{v',v,m)p{x,v,m,t)dv', 

Jv 

(3.8) 


poo p 

— AS + S = h(x, t) := / / p{x,v,m,t)dvdm, 

Jo Jv 

(3.9) 

and the initial data 

p = 0 . 

(3.10) 

From ()3.9p we write S 

= G * h. By Young’s inequality. 




(3.11) 


Next, from (13.81) we have 

dtp + v ■ VxP + F{m, 8 i)dmP 

= -dmF{m, 8 i)p-dm(p 2 j dsF{m,{l - e )82 + 98i)d6'js 
+ / T{v,v',m)p{x,v',m,t) — T{v',v,m)p{x,v,m,t)dv'. 

Jv 

Integrate along the characteristics dehned by (j2.10p from 0 to t, with 8 replaced by Si, we get 

p{x,v,m,t) =po(X(0),u,M(0)) - f 5mF(M(r), Si(X(t), r))p(X(T), u, M(t), r)dr 

Jo 

- dm{p 2 1^' dsF{m, (1 - 6)82 + 98i)d9) (X(r), u, M(r), t)S(X(t), r)dr 

+ f [ r(u, u', M(r))p(X(T), u', M(r), r) — r(u^, u, M(r))p(X(T), u, M(t), r)(iu^(ir. 

Jo Jv 
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Take ^ ^ norm on both sides, we obtain 

l|p(i)||Li_„^^ < ||po||li,„,^ +[n(maxJ5(X(s),s)|) + 2FdCr] / 

+ (||pIIl-(lJ.^) + |K||Ljc(ii^^))n(maxJ5(X(s),s)|) ||5||ii. 

Notice from Proposition 13.31 that p,Pm £ (I3.1ip . we get 

IIpWIIli,„,^ < IIpoIIli,„,^ + C(t) f \\p{s)\\li ds, 

Jo 


where C{t) is bounded for any given t. Recall the zero initial data (I3.10I1 . thus by applying Gronwall’s 
inequality we have 

IIpWIIli _ = 0 , 


then we have proved the uniqueness of the global weak solution. 


4 Fast adaptation limit 

In this section we investigate, as e ^ 0, the limiting behavior of the system 


□ 


dtp^ + V ■ VxP^ + dm[ 


F{m,S^) 

e P\ = 


V 


T{v, v', m)p^(x, v', m, t) — T{v', v, m)p^{x, v, m, t)dv' 


roo r 

'{x,t) := / / p’^{x,v,m,t)dvdm, 

Jo Jv 


-/\S^ + S‘^ = n^( 

>0 Jv 

p^(x, u, m = 0, t) = 0, p’^ix, v,m = +oo, t) = 0, 

(4.1) 

with initial data po which satisfies (ll.2p - (|1.3l) . Below we denote A4(ll) the space of Radon measures 
on fl, Co is the Banach space of continuous functions which vanishe at oo, and the notation 


tm. 


roo 

p^{x,v,t) := / p^{x,v,m,t)dr 

Jo 

First of all, we recall the a priori estimates including the parameter e, which is 
Proposition 4.1 Under the same assumptions in Theorem \2.21 The solution to m satisfies 

l|P^('i ■)'!^)llL°°(R^xyxR+) — IIpo||l°° (l + —e E ), Vt > 0, 

IIp^(') "j ■) ^)llLi(RdxyxR+) ~ IIpoIIlIj Vt > 0, 

\\S^i;t)\\L^nL-iRi) < \\Po\\L^+Vfi\po\\L^il + 2VdCTte^^‘^^rt^, Vt > 0, 
\\n^;t)\\L^iRi) ^ Vd\\po\\L^{l + 2VdCTte^^'^^rf^, Vt > 0, 
\\P%-vU)\\L°-{RixV) < ||po||l-( 1 + 2RrfCrte2^''‘^'^*), vt > 0. 


and 


(4.2) 

(4.3) 

(4.4) 

(4.5) 

(4.6) 
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From the above proposition, we have 


Lemma 4.2 Let he solution of (j4.ip and further assume (ll.6p - (ll.7jl . Fix any T > 0. Then, 


up to a subsequence, we have 

p'^ ^ p in X y X X [0, T]) — u;*, (4.7) 

p^^p in L°°{RixV x[0,T])-w*, (4.8) 

poo 

p{x,v,t) = / p{x,v,m,t)dm, (4.9) 

Jo 

S'^{x,t) ^ S{x,t) strong in L'^(M^ x [ 0 ,T]), VI < 5 < oo, (4-10) 

and further 

S%x,t) ^ S{x,t) in Co{Rix[0,T]). (4.11) 


We postpone the proof of Lemma 14.21 to the end of this section. The main theorem in this section is 

Theorem 4.3 Under the assumptions in Lemma \4-^ Let (p, p, S) he defined in Lemma \4-^ Then 

p{x, V, m, t) = p{x, V, t)5{m — mo(5)), (4.12) 

where mo(5') is defined in (|1.6I) such that F{mo{S),S) = 0. And {p,S) satisfies, in the sense of 
distribution, 

dtp + v-VxP= f T{v,v',moiS))p{x,v',t)-T{v',v,mo{S))p{x,v,t)dv', (4.13) 

Jv 

— AS + S = n{x,t)-.= p{x,v,t)dv. (4.14) 

Jv 

The initial data of the resulting limit equation (I4.13p is given by pQ = pQ{x,v,m)dm. 

Proof. Multiply Equation (I4.1h i by e to get 

dm[F{m, S^)p^] = £ [T{v,v',m)p^{x,v',m,t) — T{v',v,m)p^{x,v,m,t)]dv' — sdtp^ — ev ■ VxP^■ 

Jv 

Recall we had ()4.7p that p^ ^ p weak* in A4(]R^ x R x M+ x [0, T]), and ()4.1ip that S^{x, t) S{x, t) 
in C'o(M^ X [0, T]), therefore we can pass to the limit in weak sense as e ^ 0, up to a subsequence, to 
get 

dm[F{m, 5')p] =0 in the sense of measure, (4-15) 

thus 

F{m, S)p = constant = 0 in the sense of measure. 

Using the assumption ()1.6p on F, also using the fact ()4.9p . we readily conclude (I4.12p . 

On the other hand, take the integration of (l4.1D i with respect to ni we have 

dtp^ + V ■ VxP^ = / / T{v,v',m)p^{x,v',m,t) — T{v',v,m)p^{x,v,m,t)]dv'dm. 

Jo Jv 


14 









By (14.81) . we pass to the limit as e ^ 0 up to a subsequence, and find, in the sense of distribution, 

/*oo p 

dtp + v-'S/xp= / / T{v,v',m)p{x,v',m,t) — T{v',v,m)p{x,v,m,t)]dv'dm. 

Jo Jv 

Recall that p has the form in (|4.12l) , then the equation of p is thus (14.131) . 

Next, (I4.14j) is derived by passing to the limit as e —>• 0 up to a subsequence on (14.11) 9. using the 
property (14.101) . □ 

Proof. Proof of Lemma m 

1. Recall the a priori estimate, {p^} is bounded in x R x M+ x [0,T]), thus w* precompact 

in A4(M^ X R X M+ x [0, T]), so we have (14.7p . 

2. Similarly, recall (14.61) . {p'^} is uniformly (with respect to e) bounded in x R x [0,T]) thus 

w* precompact so we have ()4.8h . 

3. To prove ()4.9p . we have to control the tail for large m. We notice from the assumption ()1.6p that 

F{m,S^) < 0 when m > m+, 

where m+ is a uniform upper bound of mo{S^). Now we take a smooth non-decreasing function (p{m) 
such that 


(p{m) = 0 when m < m_|_, p{m) < m when m+ < m < 2m_|_, p{m) = m when m > 2m_|_, 
multiply with the equation (|4.1h i . and integrate with respect to x,v,m, then we have 

F{m,S^ 


thus 


"W / [ f p{'m)p^dvdxdm = [ [ [ ^ p^dvdxdm < 0, 

dtJm+kiJv Jm+kiJv ^ e ^ 

roQ r p 1*00 p p poo p p 

/ / / mp^dvdxdm < / / p{m)p‘^dvdxdm < / / ip{m)podvdxdm 

J2m+ J^i Jv Jm+ Jm.^ Jv Jm+ Js.^ Jv 


< 


poo p 

Jo Jm 


'v 


mpodvdxdm < oo. 


(4.16) 


The last inequality is based on the assumption (II.7p and gives the control at oo. Now choose a test 
function (p{x, v, t) G Cq and a smooth cutoff function 


0 < Xr{i^) ^ = 1 for m < R, Xr{ir) =0 for m > 2R. 


Then we compute 

rT 


/ / / 4i{x,v,t)p^dvdxdt 

Jo JmiJv 



4>{x, V, t)p^dmdvdxdt 
(j){x, V, t)xR{'m)p^dmdvdxdt 

OO 

(j){x, V, f)(l — XR{iR))p^dmdvdxdt. 


(4.17) 


15 












For the first term on the right hand side, note from ()4.7I) that ^ p weak*, for any given R, 


pT p p poo pT p p poo 

/ / / / (j>{x,v,t)xR{'m)p^dmdvdxdt ^ / / / / 4>{x,v,t)xR{'m)pdmdvdxdt 

Jo Jmi Jv Jo Jo Jm Jv Jo 


as E —y 0. 


For the second term, the control at oo estimate (I4.16P ensures 


/O J'R.i.JVJO 


p poo pT p p poo 

/ / (p{x,v,t){l — XR{'m^))p‘^dmdvdxdt < \\4>\\co / / / / p^dmdvdxdt 

Jv Jo Jo Js-i Jv Jr 


<ii 


I Co 


If 


R Jo JRi Jv Jr 

thus let i? —)• oo in (j4.17p we have 

cT r r rT 


mp^dmdvdxdt —)• 0 as i? —>• oo. 


fJ r r f f 

/ / / (j){x,v,t)p‘^dvdxdt ^ / / / / 4>{X}V,t)pdmdvdxdt 

Jo JRi Jv Jo Jm.i Jv Jo 


as E —y 0. 


Since we already have (14.81) . by uniqueness of the limit we proved (|4.9p in the sense of measure. 

4. For {S''^(x,t)}, recall (|4.4I) . it is bounded in x [0,T]) uniformly with respect to e, to 

show strong compactness, we use the claims 


Local compactness in space: 

e Vte[0,T], Vl<g<oo, 

Local compactness in time: 

Vte[0,r], Vl < g < oo. 

Control at oo: 


I {x){S‘^ydx < oo, Vt G [0,T], VI < g < oo. 
jRi 


(4.18) 


(4.19) 


(4.20) 


The first two claims are elliptic estimates which can be similarly derived as in last section. We only 
prove P4.20I) . Recall the equation for p^: 

dtp’^ + V ■ + dm[ -^- P^] = / T{vi v', m)p‘^{x, v', m, t) — T{v', v, m)p’^{x, v, m, t)dv'. 

e J V 

Multiply both sides of the equation by (x) and integrate with respect to x, v, m, we have 


jL 

dt 


[ r [ {x)p^ 

jRi Jo Jv 


dvdmdx = 


V ■ X 


p^dvdmdx < IIpoIIlc 


that is, 


thus we obtain 


yo Jv (x) 

^ j (x)n^dx < IIpoIIlL 
dt jRi 

/ {x)rfdx is uniformly bounded for all t G [0,r]. 

jRi 


(4.21) 
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Next, recall the equation for S^: 


-AS^ + S^ = 

multiply both sides of the equation by , for any q large and fixed, then integrate with respect 

to X, we have 

[ {x){S^y-H-^S^ + S^)dx = [ {x){S^f-^rfdx, 

J'Ri Jm.i 


then 


Note here that 


[ {x){S^)^dx+ [ V{{x){S^)^-^)-VS^dx= [ {x){S^y-^n^dx. 
jRi jRi JRi 


[ V{{x){Sy‘^-y-VS^dx = [ {x){p-l){Sy'^-^VS^-VS^ + [ • V(x) 

jRi Jr^ Jr^ 

= [ (x)(p-l)(5^r2|V5^|2- / -(5 TA(x) 

jRd jRd q 

= f w(p-i)(s'r"ivs'p-/ V)»(^+7^). 

jRi JRi q {x) {xy 

Then, we find 

[ {x)isy'^dx+[ {x){q-i)isy'^-^\vs^\^dx= [ -(syy^^ + -^)dx+ [ (x)(5^)''-Vdx, 

JRi JRi JRi q \x} {x}'^ jRi 


that is 


[ {x){S^ydx < - [ (5^)^dx + ||S''^|| [ {x)n^dx, 

JRi q JRi "" JRi 


using again (I4.4h that {5''^(-, t)} is uniformly bounded in L^. n for all t € [0, T], and the fact (I4.21jl . 
we conclude, for any fixed q large, 



< oo, Vt e [0,T], 


Note also ()4.20l) already holds for q = 1 , which is because we can use the similar argument as for 
(|2.16p . therefore we have proved (14.2011 for all 1 < g < oo. 

In conclusion, by using again the strong compactness criterion in Brezis [8], we have (j4.inp from the 
above claims. 

5. To prove ()4.1ip . we first estimate the control of at x = oo. For simplicity of notation, we 
abbreviate the time variable in the following. 

Fix R > 1 . For |x| > R, we compute 


^|i/d^£(x) _ J \x\^^'^G{x — y)n'^{y)dy 

< j \x - y\^/‘^G{x - y)nyy)dy + j \y\^/‘^G{x - y)nyy)dy 
= : F + IF. 
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For P, note that G is only singular at the origin, we cutoff the singularity and estimate 


f \x - y\^^'^G{x - y)n%y)dy + f \x - y\^^'^G{x - y)n%y)dy 
J\x—v\<l J\x—v\>l 


F = 

/ 

\x-y 


J\x-y\<l 


< 

f 

G{x - 


J\x-y\<l 


< 

TOO 

M M -^-'x 



l^l>i 


'\x-v\>l 


\x-y\>l 

X-y\^ldG‘^{x- 


P IIl2 


recall (j4.5p . note also the bound of {n^} is trivial from conservation of mass, thus {n^} is uniformly 
bounded in L^. n L!^, and note the property (P3) of Bessel potential G{x) in Proposition 12.11 the 
integration above is finite. Then we conclude that P is uniformly bounded (the bound is independent 
of e). 

For IP, we use Holder’s inequality, with q = and q' = d, to get 


IF = 


< 


\y\^/'^G{x - y)rf{y)dy < ||G|| 


lx=\ 


||G||^^(||n=||^J J \y\n^y)dyy. 


\y\n^y) 



(4.22) 


Here we use the property of Bessel potential (P2) in Proposition 12.11 together with the estimates (I4.5|l 
and ()4.2ip . we conclude from (j4.22l) that IF is also uniformly bounded (the bound is independent of 
e). Combine the above, we conclude 

I^|l/d^e(^) ^ ^ 1^1 y Jl (4.23) 


Next, by Morrey type Sobolev embedding Theorem (see Evans m), note the elliptic estimates 
(j4.18p - ()4.19p imply that the convergence of S^{x,t) S{x,t) is actually in Holder space x 

[0,T]), for any 6 G (0,1). Together with the estimates (|4.23l) . which yields that S{x) vanishes at 
infinity, then (14.lip is proved. □ 


5 Conclusion 

We have considered a nonlinear kinetic chemotaxis model with internal dynamics incorporating signal 
transduction and adaptation. Under some quite general assumptions on the model and the initial 
data, we have proved the global existence of weak solution by using the Schauder fixed point theorem. 
More precisely, for our mathematical treatment, we generalise the assumptions on the adaptation rate 
F and turning kernel T in the model than that in [MSI IS]. Compare with the global existence 
result in |14] . our result holds for any physical space dimensions. Moreover, the uniqueness of weak 
solution is also derived, based on some further regularity estimates on the solutions, following the 
method devised in [26] . 

Next, we considered the fast adaptation limit of this model to Othmer-Dunbar-Alt type kinetic 
chemotaxis model. This limit gives some insight to the molecular origin of the chemotaxis behaviour, 
by incorporating information about microscopic intracellular processes such as signal transduction and 
response into the chemotaxis description. This was done in Hade] for a highly simplified description of 
intracellular dynamics, where linear dynamics for the response to an extracellular signal was assumed. 
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We remark also that in order to derive the molecular origin of the chemotaxis behaviour, we did not 
use any moment closure to derive the closed evolution equation for the macroscopic density of cells as 
in [laiis], instead, a kinetic type limit equation ()4.13p is arrived, with turning kernel incorporating 
information from microscopic intracellular processes. 

In our analysis, the fast adaptation limit is derived by extracting a weak convergence subsequence 
in measure space. For this limit, the first difficulty is to show the concentration effect on the internal 
state. When the small parameter e, the adaptation time scale, goes to zero, we prove that the solution 
converges to a Dirac mass in the internal state variable. Another difficulty is the strong compactness 
argument on the chemical potential, which is essential for passing the nonlinear kinetic equation to 
the weak limit. 

For future works, it is interesting to consider the case with several chemical reactions as related to 
the model proposed in Erban-Othmer [mns]. We also mention that other types of scaling and related 
limits can be considered. For example, based on both fast adaptation and stiff response, the paper m 
studied how the path-wise gradient of chemotactic signal arises from intra-cellular molecular content. 
Several other rescallings are possible. For instance, an open problem is the asymptotic behaviour of 
the hyperbolic scaling, with (jl.lh i replaced by 


Firn 1 f 

dtp + V ■ VxP + dm[ -^— p] = - T{v, v', m)p{x, v', m, t) — T{v', v, m)p{x, v, m, t)dv'. 

£ ^ Jv 


The hyperbolic limit is quite different from the fast adaptation limit, because, in that case, even the 
a priori estimate for p is not uniform in e, so the limit equation is unclear, a difficult part of that will 
be to determine how the extracellular signal feeds into the chemotaxis response of the cells. 
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